Selectivity of the Nucleon Induced Deuteron Breakup and Reiativistic Effects 



H. Witala, J. Golak, and R. Skibinski 

Institute of Physics, J agiellonian University, PL-30059 Cracow, Poland 
(Dated: February 9, 2008) 

Theoretical predictions for the nucleon induced deuteron breakup process based on solutions of 
the three-nucleon Faddeev equation including such reiativistic features as the reiativistic kinematics 
and boost effects are presented. Large changes of the breakup cross section in some complete 
configurations are found at higher energies. The predicted reiativistic effects, which are mostly of 
dynamical origin, seem to be supported by existing data. 

PACS numbers: 21.30.-x, 21.45.+V, 24.10.-i, 24.70. +s 



in ■ 

o : 
o ■ 

<N ■ 

Cu- 
te ■ 

m ; 
o ■ 



> 
oo 
On 
O 
<3\ 
O 

in 
o 



o 

=5 



X 



Recent studies of elastic nucleon-deuteron (Nd) scat- 
tering and nucleon induced deuteron breakup revealed 
a number of cases where the nonrelativistic description 
based only on pairwise nucleon-nucleon (NN) forces is 
insufficient to explain the three-nucleon (3N) data. This 
happens in spite of the fact, that these high precision NN 
potentials describe very well the NN data set to about 350 
MeV laboratory energy. Those findings extended explo- 
ration of the properties of three-nucleon forces (3NFs) to 
the reactions in the 3N continuum. Such forces appear 
for the first time in the 3N system where they provide 
an additional contribution to a predominantly pairwise 
potential energy of three nucleons. Generally speaking, 
the studied discrepancies between a theory based only 
on NN potentials and experiment become larger for in- 
creasing energy of the 3N system. Adding a 3N force to 
the pairwise interactions leads in some cases to a better 
description of the data. The best studied example is the 
discrepancy for the elastic angular distribution in the re- 
gion of its minimum and at backward angles 0, 0, Q 
This clear discrepancy can be removed at energies be- 
low w 100 MeV by adding modern 3NFs to the nuclear 
Hamiltonian. Such a 3NF, mostly of the 27r-exchange 
character, must be adjusted individually with each NN 
potential to the experimental binding energy of 3 H and 
3 He. At energies higher than ks 100 MeV current 3NFs 
improve only partially the description of cross section 
data and the remaining discrepancies, which increase 
with energy, indicate the possibility of reiativistic effects. 
The need for a reiativistic description of 3N scattering 
was also raised when precise measurements of the to- 
tal cross section for neutron-deuteron (nd) interaction Q 
were analyzed within the framework of the nonrelativistic 
Faddeev calculations 5] . Also there the NN forces alone 
were insufficient to describe the data above « 100 MeV. 
The effects due to reiativistic kinematics considered in 
ref. [f| were comparable at higher energies to the effects 
due to 3NFs. This demonstrates the importance of a 
study taking reiativistic effects in the 3N continuum into 
account. 

Recently first results of reiativistic 3N Faddeev calcu- 
lations for the elastic Nd scattering have become avail- 
able 0. The reiativistic formulation applied was of 
the instant form of reiativistic dynamics 7]. A start- 
ing point of this formulation for 3N scattering is the 



Lorentz boosted NN potential V(k, k ';P ) which gener- 
ates the two-nucleon (2N) boosted t matrix t(k, k '; P ) in 
a moving frame via a standard reiativistic 2N Lippmann- 
Schwinger equation 
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t(k, k';P) = V(k, k';P) + j d 3 k" 
V(k,k";P)t(k",k';P ) 
y/(2u(k')) 2 + P 2 - \l{2to{k ")) 2 + P 2 + le 



The NN potential in an arbitrary moving frame V{P ) 
is obtained from the interaction v defined in the two- 
nucleon cm. system by |8j 



V(P) = yj(2u(k) + v) 2 + P 2 - yj(2u(k)) 2 + P 2 . (2) 

The reiativistic kinetic energy of three equal mass (m) 
nucleons in their cm. system can be expressed by the 
relative momentum k in one of the two-body subsystems 
and momentum of the third nucleon q (the total momen- 
tum of the two-body subsystem is then P = —q) as 
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k 2 is the momentum dependent 



where 2a; (k) = 2\fm? 
2N mass operator. 

The Nd scattering with neutrons and protons interact- 
ing through a NN potential V alone is described in terms 
of a breakup operator T satisfying the Faddeev-type in- 
tegral equation [!J, 



T\cj}> = tP\<f>> +tPG T\4>> 



(4) 



The permutation operator P = P12P23 + P13P23 is given 
in terms of the transpositions P^ , which interchange nu- 
cleons i and j. The incoming state \<f> >= \qo > \(pd > 
describes the free nucleon-deuteron motion with the rela- 
tive momentum q$ and the deuteron wave function \(j>d >■ 
The Go = E+i \_ Ha is the free 3N propagator with total 
3N cm. energy E expressed in terms of the initial neu- 
tron momentum qo relative to the deuteron 
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where M4 is the deuteron rest mass. 

The transition operators for elastic scattering, U, and 
breakup, Uq, are given in terms of T by 

U = PG^ 1 + PT 

U = (1+P)T. (6) 

The state Uo\(j> > is projected onto the state \4>q > which 
describes the free motion of the three outgoing nucleons 
in the 3N cm. system in terms of the relative momentum 
of the 2N subsystem A^jy c .m.(2 — 3), defined in the 3N 
cm., and momentum of the spectator nucleon q defined 
above: \(f>o >= \k 3 N cm. (2 — 3) q >i. This leads to the 
breakup transition amplitude 

i 

The choice of the relative momentum k in the NN 
cm. subsystem and the momentum q of the specta- 
tor nucleon in the 3N cm. system to describe config- 
uration of three nucleons is the most convenient in the 
relativistic case. In the nonrelativistic limit the mo- 
mentum k reduces to the standard Jacobi momentum 
p 01 • To solve Eq. (0} numerically partial wave decompo- 
sition is still required. The standard partial wave states 
\pqa >= \pq(ls)j(X^)IJ(t^)T > TGj, however, are gen- 
eralized in the relativistic case due to the choice of the 
NN-subsystem momentum k and the total spin s both 
defined in the NN cm. system. This lead to Wigner spin 
rotations when boosting to the 3N cm. system 0,0; re ~ 
suiting in a more complex form for the permutation ma- 
trix element 6] than used in the nonrelativistic case ^3 • 
A restricted relativistic calculation with j < 2 partial 
wave states showed that Wigner spin rotations have only 
negligible effects 6] . Due to this we neglected the Wigner 
rotations completely in the present study. To achieve 
converged results at energies up to « 250 MeV all par- 
tial wave states with total angular momenta of the 2N 
subsystem up to j < 5 have to be used and all total an- 
gular momenta of the 3N system up to J = 25/2 taken 
into account. This leads to a system of up to 143 coupled 
integral equations in two continuous variables for a given 
total angular momentum J and total parity ir = (— )' +A of 
the 3N system. For the details of our relativistic formula- 
tion and of the numerical performance in the relativistic 
and nonrelativistic cases we refer to Ref. 0, |j| . 

In the present study we applied as a dynamical input a 
relativistic interaction v generated from the nonrelativis- 
tic NN potential CDBonn according to the analyti- 
cal prescription of ref. |l2j |. This analytical transforma- 
tion allows to obtain an exactly on-shell equivalent to the 
CDBonn relativistic potential v which provides the cor- 
responding relativistic t matrix. The boosted potential 
was not treated in all its complexity as given in ref. [l3T ] 
but a restriction to the leading order term in a P/u and 
v/oj expansion was made 

V(k,k';P) = v(k,k') 
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FIG. 1: Differential cross section for elastic Nd scattering at 
250 MeV. The dashed and solid lines are the results of the 
relativistic and nonrelativistic Faddeev calculations based on 
the NN potential CDBonn. The pd data are from [l9||. 
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8Vm 2 + k 2 Jm 2 + (k') 2 

The quality of such an approximation has been checked 
by calculating the deuteron wave function 4>d(k) of the 
deuteron moving with momentum P for a number of 
values corresponding to incoming nucleon lab. energy 
< 250 MeV. The resulting deuteron binding energies and 
the deuteron D-state probabilities for the deuteron in mo- 
tion are close to the values for the deuteron at rest. 

In Fig. n we show the nucleon angular distribution for 
elastic nucleon-deuteron scattering at Ef ab = 250 MeV. 
It is seen that, like in a study of Nd elastic scattering 
in ref. @ where the AV18 [llj NN potential instead of 
CDBonn have been used, relativistic effects for the cross 
section are restricted only to the backward angles where 
relativity increases the nonrelativistic cross section. At 
other angles the effects are small. In spite of the fact that 
the relativistic phase-space factor increases with energy 
faster than the nonrelativistic one (at 250 MeV their ratio 
amounts to 1.175), the relativistic nuclear matrix element 
outweighs this increase and leads for the cross section in a 
wide angular range to a relatively small relativistic effect. 

The breakup reaction with three free outgoing nucle- 
ons in the final state provides a unique possibility to ac- 
cess the matrix elements of the breakup operator T with 
specific values of momenta \k\ and |g| in a pointwise man- 
ner. Each exclusive breakup configuration specified com- 
pletely by 3N cm. momenta fej of outgoing nucleons re- 
quires three matrix elements j < k(kj,kk),q — ki\T\<j) > 
with k) = (1, 2,3) and cyclical permutations, and k 
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and q providing the total 3N cm. energy 



E = dA(m 2 + k 2 ) 



q 2 + \J m 2 



q 



(9) 



This is entirely different from the elastic scattering where, 
due to continuum momentum distribution of nucleons in- 
side the deuteron a broad range of fc|- and |<j]-values 
contributes to the elastic scattering transition matrix el- 
ement. That particular selectivity of the breakup singles 
out this reaction as a tool to look for localized effects 
which when averaged are difficult to see in elastic scat- 
tering. 

This selectivity of breakup helps to reveal relativis- 
tic effects in the 3N continuum. Even at relatively low 
incoming nucleon energy Ej^ b = 65 MeV they can be 
clearly seen in cross sections of some exclusive breakup 
configurations as examplified in Figs. and [3] For the 
configuration of Fig.|2the angles of the two outgoing pro- 
tons detected in coincidence were chosen in such a way 
that for the arc-length S « 30 MeV all three nucleons 
have equal momenta which in the 3N cm. system lie in 
the plane perpendicular to the beam direction (symmet- 
rical space star (SSS) condition). For the configuration 
from Fig. [3| at the value of S » 46 MeV the third, not 
observed nucleon is at rest in lab. system (quasi-free scat- 
tering (QFS) geometry). In these two breakup configura- 
tions the inclusion of relativity lowers the cross section: 
by « 8% in the case of SSS and by w 10% in the case of 
QFS. In the lower parts of Figs. Eland [3| contributions 
to this effect due to kinematics and dynamics are shown. 
The five-fold differential cross section can be written as 



d 5 a 



= ( 



E 



<0o|W> I') Pkin, (10) 



with the kinematical factor pki n containing the phase- 
space factor and the initial flux. The transition proba- 
bility for breakup | < <f>o\Uo\<f> > | 2 , averaged over the 
initial m,i n and summed over final m ou t sets of parti- 
cles spin projections, forms the dynamical part of the 
cross section. In the lower parts of figures the ratio of 
the relativistic to the nonrelativistic kinematical factor 
Pkin I Pkin as a function of S is shown by the dashed 
line. The corresponding ratio for the dynamical parts 
of the cross section is shown by the solid line. As seen in 
Fig. for the QFS configuration the whole effect is due 
to a dynamical change of the transition matrix element. 
For this configuration the nonrelativistic and relativistic 
kinematical factors are practically equal for large region 
of S- values (see Fig.^J. For SSS about 30% of the total 
effect is due to a decrease of the relativistic kinematical 
factor with respect to the nonrelativistic one (see Fig.[5J). 

The cross sections in these particular configurations 
are rather stable with respect to exchange between mod- 
ern NN forces, combining them or not with three-nucleon 
forces Due to that relativistic effects seem to ex- 

plain the small and up to now puzzling overestimation of 
the 65 MeV SSS cros section data 16| by modern nuclear 




FIG. 2: Five-fold differential cross section d 5 a / d^lid^dS 
as a function of arc-length S for the SSS configuration (see 
text) of d(p,pp)n reaction at 65 MeV. The lab. angles of 
the outgoing protons are 9i — 62 — 54° and $12 = 120°. The 
solid and dashed lines are the results of the nonrelativistic and 
relativistic Faddeev calculations based on the NN potential 
CDBonn. The 65 MeV pd data are from 16j. The lower part 
shows the ratios of the relativistic and nonrelativistic phase- 
space factors (dashed line) and of the cross sections divided 
by the corrresponding phase-space factor (solid line). 
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FIG. 3: The same as in Fig. |5] for the QFS geometry (see 
text). The lab. angles of the outgoing protons are 81 = 62 = 
44° and <j> X2 = 180°. The 65 MeV pd data are from 0. 



forces and can account for the experimental width of this 
QFS peak 

At higher energies selectivity of breakup allows us to 
find the configurations with significantly larger relativis- 
tic effects. In Fig.Hthis is examplified at Ef ab = 200 MeV 
and the predicted effects of up to w 60%, which are 
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mostly of dynamical origin, seem to be supported by the 
data of ref. 0. 




oica ; l t I I I I I I I I I I I I I I 




: r 

40 80 80 100 

E, [MeV] 



FIG. 4: Five-fold differential cross section d 5 a / dQidfadEi 
as a function of energy of the detected proton Ei for d(p,pn)p 
reaction at 200 MeV. The polar lab. angle of the outgoing 
proton and neutron is 9\ = 52° and 82 = 45°, respectively. 
The azimuthal angle 4>i2 = 180°. For the description of lines 
and the bottom part see Fig. |2] The data are from fTs| . 

The selectivity of complete breakup is gradually lost 
when uncomplete reactions are considered. In the total 
nd breakup cross section the effects disappear. Integrat- 
ing over all available complete breakup configurations 
provides nearly equal relativistic (90.25 mb at 65 MeV 
and 43.37 mb at 250 MeV) and nonrelativistic (91.12 mb 



and 45.41 mb) total breakup cross sections. Also inte- 
grated elastic scattering angular distribution (71.25 mb 
and 9.33 mb - relativistic, and 71.40 mb and 9.57 mb 
- nonrelativistic) and the total cross section for the nd 
interaction do not reveal significant relativistic effects. 
This shows, that the discrepancies between theory and 
data found in previous studies at higher energies for the 
total cross section and elastic scattering angular distri- 
butions, which remain even after combining NN poten- 
tials with 3NFs, have to result from additional contri- 
butions to the 3N force, which have different than the 
27r-exchange character. 

Summarizing, we showed that selectivity of the com- 
plete breakup reaction enables us to reveal in 3N con- 
tinuum clear signals from relativistic effects. Existing 
breakup data seem to support the predicted effects, when 
the relativity is included in the instant form of relativistic 
dynamics as proposed by Bakamjian and Thomas. Pre- 
cise complete breakup data at energies around 200 MeV 
are welcome to further test these predictions. The QFS 
breakup configurations due to their large cross sections 
and insensitivity to the details of nuclear forces are fa- 
vored for this purpose. 
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